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Abstract. We present a general formalism to derive the evolution equations 
describing ID and isotropic 2D interface-like systems, based on symmetries, 
conservation laws, multiple scale arguments, and exploiting the relevance of coarsening 
dynamics. Our approach becomes specially significant in the presence of surface 
morphological instabilities and allow us to classify the most relevant nonlinear terms 
in the continuum description of these systems. The theory applies to systems ranging 
from eroded nanostructures to macroscopic pattern formation. In particular, we show 
the validity of the theory for novel experiments on ion plasma erosion. 
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1. Introduction 

Many surfaces and interfaces, which are very different in nature and that occur at 
very different lengthscales, produce interfaces which are starthngly similar. This 
similarity originates in the competition between stabilizing and destabilizing physical 
mechanisms, that are insensitive to the specific value of the interface height h{x, t) 
(this function providing the height of a interface above position x on a ID substrate, 
at time t). Namely, such mechanisms are invariant under global height translations 
of the form h{x,t) h{x,t) + C,, with C, an arbitrary constant (shift simmetry [1]). 
For instance, interesting [(sub)micrometric] surface features develop by growth or 
erosion using various techniques, such as electrochemical deposition (ECD) [2j, chemical 
vapor deposition (CVD) [3], ion beam sputtering (IBS) [1] or molecular beam epitaxy 
(MBE) [5], but remarkably also in macroscopic systems, such as aeolian sand dunes [6] 
or underwater (vortex) ripples in sand [7]. 

For cases in which dynamical instabilities are irrelevant or absent, the shift 
symmetry has been argued to lead, in the presence of noise, to scale invariant surface 
morphologies. This is the generic scale invariance [8] associated with the universality 
classes [9] of kinetic roughening. However, the situation differs in the presence of 
morphological instabilities, that are the landmark of pattern formation. In these cases, if 
the height field is seen as a measure of the amplitude of perturbations around a reference 
homogeneous state, the conservation law associated with the shift symmetry leads to 
a large-scale instability [in]. This might prevent a unified description by means of 
a universal equation, such as the Ginzburg-Landau equation is for the case of short- 
wavelength instabilities [TOl [H]. Thus, although the systems mentioned above are 
seemingly governed by similar equations, a general theory of evolving interfaces in which 
patterns arise is still lacking. In particular, it seems that for every experimental system 
a specific theory needs to be developed. 

In this paper, we present a general formalism for pattern formation in ID and 
isotropic 2D surfaces, that relates with previous approaches in the theory of dynamic 
scaling of rough interfaces [H [121 [13] , exploiting systematically the formulation of the 
interface evolution within the assumption of a (generalized) relaxational dynamics. 
Through symmetries and multiple scale arguments, our approach helps to understand 
why some effective equations presenting a morphological instability appear almost 
ubiquitously in the experimental systems mentioned above. Our method enables 
classification of the most relevant nonlinear terms for these systems, and stresses the 
relevance of the phenomenon of coarsening to their continuum description. We have also 
performed some experiments on ion plasma erosion to explore non-trivial implications 
of our theory. Moreover, we have applied it to other relevant and diverse fields, ranging 
from thin film production to surface patterns in macroscopic systems, such as fluid waves 
or some granular systems. 

In principle, we consider that the evolution of h is local and, mathematically, takes 
the form dth = G{{h}) + rj, where {h} stands for h and its spatial derivatives, and rj 
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is a noise term that we will neglect hereafter. Previous phenomenological proposals for 
the functional Q have been built upon geometry arguments [12], or upon symmetry and 
conservation requirements [9]. In the first case, the evolution is governed by the local 
curvature and its derivatives along the interface. This excludes contributions to Q that 
are due to, e.g., anisotropies [H]. Likewise, the second approach can be misleading. 
For instance, if one considers non-conserved interface dynamics (namely, if Q cannot 
be derived from a current) under the shift symmetry, one could naively drop terms 
like \d'^h^ (because it can be written as dxj) or hd^h (because it is not shift invariant). 
However, Q could contain a term proportional to {dxhy, namely, the sum of both terms. 
This maybe innocuous in the presence of generic scale invariance, but not in the presence 
of instabilities, hence a systematic method to classify the terms in the equation is needed. 



2. Theory 

2.1. N on- conserved dynamics 

We proceed on by studying separately systems with non-conserved and conserved 
dynamics. As in Refs. [HllS], we assume that the dynamics of the systems is relaxational 
in a generalized sense (namely, with a non-constant mobility). Although this assumption 
is not strictly necessary, it clarifies greatly the interpretation of the coefficients of the 
non-linear terms below. Hence, 

e,, . -vm/-^, (1) 

where F (> 0) and F are, respectively, a function and a functional of the height field 
and its spatial derivatives. Hereafter, we restrict {h} to be {h,dxh,dlh,dlh} |T5]. 

If the relaxation rate F is a constant, symmetry under global shifts forbids any 
explicit dependence of F on h, otherwise, what needs to be independent of h is the right 
hand side of Eq. ([T]). Consequently, one can write (note our difference in sign convention 
from [1]): 

^(/i, dxh, d% dlh) = e-^'^J^.idxh, d% ^,h), (2) 
T{h,dxh,dlh,d^,h) =e'''V^{dxh,d%dlh), (3) 

where the height scale s is the generator of the group of symmetry under translations 
in and we have written F{{h}) = J J-'{{h}) dx. 

Close to the instability threshold, we can assume that local slopes are small, 
\dxh\ -C 1, and we can expand Fi and J^i in power series. Therefore, 

Ti = Y.^^''''\dM{dlhy{dlh)\ (4) 

-^1= f'^'''\dxh)\dlhy {dlhf . (5) 

Using these expansions in Eq. ([T]), we obtain a list of terms compatible with the shift 
symmetry. The resulting equations are quite involved. Hence, we first consider the 
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linear terms. Later we will include the most relevant non-linear terms. Thus, we find 
d^h = 14 + V\\d^h + vdlh + ^&lh - Kd^h + . . . , (6) 

where 

^± = ^7/, ^| = ^7^'°°V, (7) 
^ =s7(°^o)/ + 2s7/(°^°) + 27/('°°), (8) 
(3 =57^°°^)/, ir = -27/(°20)+27/(i°i), (9) 
and 7 = 7^''°°), / = y(ooo)^ rjj^g j^gg^j pg^j,|. ^j^g dispersion relation for single mode 
solutions of Eq. ([6]), h{x,t) ~ exp[co'gt + i/cx], reads 3fJe(co'g) = — z/g^ — Kq'^; consequently, 
a morphological instability arises if z/ < and > 0. The characteristic length- 
scale for the ensuing pattern will be roughly given by the wavelength of the mode 
maximizing Uq, namely, / = 27i^j2K/\h'\. When nonlinear terms are incorporated 
into Eq. IQ, this length-scale can grow in a coarsening process [16], that will help 
us to identify which nonlinear contributions are actually playing a role in the full 
dynamics. Close to the instability threshold z/ = 0, Eq. ([8]) can be written as 
u = —2'jf^'^^^\s — Sc)/sc = —2'yf^'^^^^e, with e a small dimensionless parameter, and 
= -27/(20o)/(^^(oio)j_^27/(oio))_ Note the length-scale / diverges as e'^/^ in the limit 
e — s> 0, as is s-independent, according to Eq. ([9]). Consequently, we can rescale length, 
time and height as x — *■ x/l, t — > t/l^ and h —>■ h/l'^, and, using the multiple scales 
method, look for the most relevant nonlinear term. Similarly to the theory of dynamic 
scaling |9], we obtain it to be A(s, s^){d^hf, with A(s, s^) = 5/7(200) _s2^j(oio)_g^j(2oo)^ 
where (s, s^) is used to stress the polynomial dependence of A on s. For many 
experimental situations, however, the hydrodynamic limit may lie beyond practical reach 
p!71 [18] , in which case the final structure depends crucially on subdominant terms, that 
must be included in the evolution equation in order to capture essential features of the 
original nonlinear system. This can be done within our multiple scale framework. For 
instance, the ratio between {dxhY and {dxh){d1h) scales as: 

{dxh){dlh) ^ 

so, in the limit / — 00, the first one dominates the second one. This calculation be 
reproduced for every term resulting in the series expansion, with the results that follow 
for the first correcting terms, in order of relevance: 

Aiio(l, s, s2)((9^./i)((9», A3oo(l, s, s'^){d^hf, 

X^{l,s)dl{dxh)\ X,o,is)idxh)idlh), (10) 
A2io(l, s, s'){d,h)\dlh), A4oo(s, s'){dxh)\ 



where 



Ano = s/7(i'i'°) + 2s7(i'0'0) /(O'l'O) + 27(1'°'°) + 3s7/(i'^'°) 
+ 3,2^/(0,0,1) + 67/(3A0)^ 

-^300 = 



S(_S^(1,0,0) f {0,1,0) _ ^(1,0,0) f (2,0,0) _ _g2^y(0,0,l) _ _g^/(l,l,0) 
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(0,1,0) j(2,0,0) 



A2,o = - 2s2[/(o.i-o)^(o,i,o) _ 3^.(0,0,1)^(1,0,0) ^ ^(2/(0'2'0) - 6/(1'°'^))] 

_^ 4^(2,0,0)^(2,0,0) _^ s(5^(l-0,0) j(l,l,0) _ 2^(0,1,0) j(2,0,0) 

+ 4/(0,1,0)^(2,0,0) ^ 3^/(2,1,0) + 2/7^2'^'°)) + 127(^'°'°) /(^'"'O^ 

+ 247/(^'0'°), (14) 



(1,0,0) /(1,1,0) _ 2^(1,0,0) /(3,0,0) _ ^2^j(l,0,l) _ _g^y(2,l,0) 




and we have used that [plhY = [dl^dxh)"^ — dxhd^h]/2. 

An example is provided by the so-called dissipation-modified Korteweg-de Vries 
equation, 



proposed [19] as the generic amplitude equation for weakly non-linear waves, appearing 
in contexts from Marangoni or Rayleigh convection to periodic waves in binary fiuid 
mixtures (see Refs. in |19j). In our surface context, it describes e.g. dynamics of (ID) 
ion-sputtered interfaces under oblique incidence, see Refs. in ^0\. In general, while the 
term {d^hY dominates asymptotically, it interrupts coarsening [211 [22]. This is actually 
the case for Eq. (fT6l) . Consequently, eventual coarsening of the pattern must be due to 
other terms in Eqs. ffTOj) . E.g., if one adds the A2 term in Eq. ffTOj) to the right hand side 
of Eq. f|T6|) . (interrupted) coarsening indeed occurs [23]. Further details and examples 
are provided below in Section [3TT1 

2.2. Conserved dynamics 

In the case of conserved dynamics for the interface, we can write [1] 



The assumption of (generalized) relaxational dynamics is now stronger than in the non- 
conserved one, the operator between F and 6F/6h constraining the form of the possible 
terms in a multiple scale expansion (see below). In order to make the problem as general 
as possible, we need to include the non- relaxational current, jnr, into the right hand side 
of Eq. ( JTTj) . The origin of this type of current has been widely discussed e.g. in the 
context of mound formation in homoepitaxial growth [2U [5] . Assuming jm. to depend 
only on derivatives of h in the form 



dth = vdlh + (3dlh - Kdlh + X{d^hf + Xuod^hd^h, 



(16) 




(17) 



jnr = T.J'^'"'\9.hy{dlhy{dlh)'', 



(18) 
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and using Eqs. dl])-® as above, Eq. (fT7|) takes the form 

dth = u'd^h + I3'dlh - K^h + X'i,oid^h){dlh) 

+ xidlid^hf + \i,,{d,hf{dlh), (19) 

where all the coefficients in Eq. (fTOll depend on f^'^^^ and J'^'^^\ Thus, 

=5%/-^^''°'°^ (20) 

f3- = - j(O'i.o) (21) 

^ jmi) + + 7 (22) 

A^.o = 25^(1'°'°)/ - J(2.0'0) (23) 

= s(s7(°'i'0)/ + 4s7/(°'i'°) + 27/(2A0))/2 (24) 

A210 = ^^(-3.7/(°'^'°) - 37/(^'°'°V2 + 37(^'°'°)//2) (25) 
An interesting result derived from the above equations is that, whenever jnr = 0, the 



condition for unstable growth reads e oc so that power counting (with I ~ e^^^^, t ^ 
and h l'^) provides terms which give a divergent scaling of the local slope, namely, 
\dxh\ scales with a positive power of /. This means that the small slope approximation 
breaks down, and the expansions above are no longer valid, in analogy e.g. with the 
so-called superrough scaling in the context of surface kinetic roughening [9], occurring 
specifically in the case of conserved dynamics. In these conditions, detailed knowledge 
of system specifics is needed and a strongly nonlinear analysis must be performed in 
order to obtain the correct interface equation, as in [25]. 

3. Results 

Thus far, we have derived the most general and relevant interface equations compatible 
with the shift symmetry. We proceed on by testing the validity of the theory and, 
more importantly, its capability to identify the right terms of the evolution equation 
in each context, based on experimental or theoretical evidences, through our procedure 
that accounts for conservation, symmetries and coarsening, in systems where a surface 
instability is present. To illustrate this method of analysis, we discuss several ID 
examples of systems with lengthscales ranging from microscopic to macroscopic, for each 
type of dynamics, both non-conserved and conserved. Moreover, some of the examples 
correspon to isotropic 2D surfaces. 

3.1. Application to one- dimensional systems 

To begin with, we consider non-conserved growth of nanometric sized patterns by Ion 
Beam Sputtering (IBS). In this technique, the surface of a solid target is bombarded 
isotropically with a flux of energetic ions. Experimentally, a typical length scale develops 
at early times, that is seen to grow later (coarsening). The system is not conserved due 
to the mass loss as a consequence of the bombarding process. Finally, the experimentally 
obtained surfaces do not develop steep slopes (so a small slope approximation is well 
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justified). With these ingredients in mind, our theory predicts that the evolution 
equation for these systems must contain the term proportional to A due to the non- 
conserved nature of the experiment; and terms A2, Aioi and A210 due to the presence of 
coarsening. Hence, 

dth = u{dlh) - K{dth) + \{d^hf + \2dl{dlh) + Aioi(9,/i)(9^/i) 

+ \2i^{d,h)\dlh) (26) 

Consistently with the above expectation, a hydrodynamic theory of erosion has been 
reported [21] which couples a fast field describing the density of mobile particles at 
the surface, with a slow one which corresponds to the observed experimental surface. 
This theory allows to relate experimental magnitudes with corresponding equation 
parameters and leads to, precisely, Eq. fl26l) . 

At intermediate length scales (microns to millimeters), we can take growth by 
ECD [2], where the surface of an aggregate grows by incorporation of cations from a 
solution. The assumption of locality above is reasonable for a small sticking probability, 
i.e., if the cations do not attach to the first visited surface site but are, rather, allowed to 
diffuse further until aggregation occurs. The cation flux at any surface point is isotropic, 
hence the system is reflection (x — —x) symmetric. Finally, the characteristic length 
of the unstable structures observed experimentally does not grow in time, i.e., there 
is no coarsening. These data suggest that surface dynamics for ECD is described by 
the Kuramoto-Sivashinsky (KS) equation, namely, V|| = /? = 0, due to the reflection 
symmetry, and that all the terms in Eqs. ( fTOl) are negligible, compatible with the absence 
of coarsening. The same scenario is observed in surface growth by CVD [3]. Indeed, 
the KS equation has been derived explicitly for both ECD and CVD from constitutive 
equations, in agreement with our results [T7] . 

Our final example of non-conserved interface growth is taken from surface 
instabilities in fluids due to surface tension gradients (see [26| for a comprehensive 
reference), an standard example being provided e.g. by Marangini-Benard convection 
|27j . In this system, a thin fluid layer is subject to an external transverse temperature 
gradient, there being non-conserved fluxes that implement the dissipation in the system. 
When heating from below a convective instability sets in whose spatial characteristics 
do not coarsen in time. This instability leads to travelling waves on the free surface of 
the system that break the x —x reflection symmetry. An asymptotic study of the 
corresponding hydrodynamic problem leads to the equation [271 ES] 

d^h = v{dlh) - K{dth) + \{d,hf + \no{d.h){dlh), (27) 

in agreement with the general expectation from our approach under the constraints 
considered. 

As an important example of conserved interface dynamics at nanometric scales, we 
consider the growth by Molecular Beam Epitaxy (MBE) of surfaces that are vicinal to 
a high symmetry crystalline orientation [28]. In this type of experiments, the molecules 
incorporated from the chamber arrive isotropically at the surface (and hence reflection 
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symmetry is fullfilled) , the relevant one- dimensional interface being given by the location 
of the step separating one terrace from an adjacent one. Under conditions in which 
the rate of adatom evaporation from the terrace is vanishingly small, the dynamics is 
conserved. It is precisely in this case that the small slope approximation has been seen 
to break down [25], as for Eq. (ITTIl for jnr = 0. However, if adatom desorption occurs on 
the terraces, the system becomes non-conserved and the relevant interfacial description 
is provided by the Kuramoto-Sivashinsky equation ^29j, as expected from our above 
arguments. 

Finally, an example of conserved interface growth in a macroscopic system is 
provided by the formation of macroscopic ripples in aeolian sand dunes [12]. In this 
case the wind direction breaks the reflection symmetry and, surface dynamics at the 
sand bed being conserved due to the conservation in the number of sand grains, one 
would expect Eq. ( |T9l) to hold as the relevant continuum equation. Indeed, this has been 
shown to be the case within the so-called hydrodynamic approach to pattern formation 
in these systems |12j . 

The results in this section, as well as some other experimental systems that can be 
analyzed along similar lines, are summarized in Table [Tj where we provide references 
to the experimental observations and to theoretical derivations compatible with our 
conclusions. We provide the reader with broader references when connection with 
specific experiments is less unambiguous. 



System 


Type 


RS 


Coars. 


Nonlinearities 


Exp. 


Theo. 


ECD/CVD 


NC 


Yes 


No 


A 


[34] 


[17] 


IBS 


NC 


Yes 


Yes 


A, A2,101,210 


[35] 


m 


Thin films 


NC 


Yes 


Yes 


A, A2 


[36] 


m 


Fluid waves 


NC 


No 


No 


A, Alio 


[26] 


[26] 


MBE (singular) 


C 


Yes 


Yes 


^210 


^ 








MBE (vicinal) 


c 


Yes 


No 


SN 


[28] 


[25] 


Sand ripples 


c 


No 


Yes 


^2.110,210 


m 


n 


Vortex ripples 


c 


Yes 


Yes 


SN 


[38J 


[7. 





Table 1. Physical examples. Non-conserved (NC), conserved (C), strongly- nonlinear 
(SN). RS means reflection symmetry {x — * —x), and "Coars." is coarsening. 

3.2. Application to novel experiments 

The examples in last subsection show how, based on general considerations, one 
can construct a general equation in good agreement with current specific theories. 
Notwithstanding, we want to emphasize the wide range of generality of our method 
by applying it to novel experiments on ion plasma erosion. Thus, a Si(lOO) wafer (2 
inches diameter) was immersed in an argon plasma, which was confined magnetically, at 
a pressure of 5 x 10~^ mbar. The only parameter that was changed in the experiments 
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was the immersion (i.e. sputtering) time. After the sputtering process, the central part 
of the wafer was analyzed by Atomic Force Microscopy operating in tapping mode with 
a silicon cantilever. Under these conditions the flux of incoming ions on the central 
part of the wafer is mainly isotropic. We plot two snapshots of the surface taken at 3h 
(Fig.[2K) and 6h (Fig. [2b). Moreover, we have measured a non-zero coarsening exponent 
n = 0.53 ± 0.02 in £ ~ t", i being the mean cell diameter. The system being isotropic, 
we can generalize our results above to this case. Out of the most relevant subdominant 
terms in Eq. ( ITOj) (once generalized to two dimensions) compatible with isotropy, A2 
[30t [22] and A210 [HI [32] are known in the literature to induce coarsening behavior when 
taken each of them in combination with Aioo and a KS dispersion relation. Extensive 
numerical simulations |23] show that the term Aioi also induces coarsening, but the term 
A400 does not. Hence, thus far we are left with an equation of the following form: 

dth = - uV^h - KV^h + X{Vhf + X2V\Vhf 

+ AioiV/i ■ V(V2/i) + A2io(V/i)2(V2/i). (28) 

^From the experimental data we can extract some quantitative information about the 
coefficients in Eq. (!28l) . Thus, a closer inspection of the experimental surface cross 
section (Fig. [1]) reveals that cell shape can be approximated by 

{uyKX,,oY^' log(cosh[(2ir/z.) V^a;]), (29) 

which is a solution of Eq. (!28l) in ID with A2 = Aioi = (and, consequently, 
approximately valid for large cells in the isotropic 2D case). This suggests that A2 and 
Aioi are negligible since large values for them would lead, rather, to parabohc cells [21] . 
Consequently, Eq. f[28l) reduces effectively to 

dth = -uV^h - KV^h + XiWhf + AioiV/i ■ V(V^/i), (30) 

the so-called convective Cahn-Hilliard equation [32]. For short times (linear regime), the 
typical lengthscale allows determination of the K/v value, while remaining parameters 
are estimated by trial and error in order to improve the resemblance with the 
experimental morphologies. Results of such numerical simulations are compared with 
experiments in Fig. [2] showing an excellent agreement between them. Thus, the small 
panels (upper insets) display two snapshots of the numerical integration of Eq. ( |30ll and 
the large panels the equivalent results from experiments. 

4. Discussion and conclusions 

Although table I is not meant to be exhaustive, it leads to some general observations. 
Thus, as anticipated above, the occurrence of coarsening is an important property in 
guiding the identification of the appropriate nonlinear terms of a given system. This 
phenomenon was associated with the seeming inability to describe patterns that are 
stable both in wavelength and amplitude by means of local height equations [TJ [16] . 
However, there may exist now counterexamples [2n [22] for this shortcoming, at least 
for an order range that is much larger than the single cell size. As a consequence, the 
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X X 

Figure 1. a) Cross section of Fig. [2^ morpliologies, both numerical (lower curve) 
obtained with v = 1, K = 0.11, A210 = 0.052 and A = —0.05 and experimental 
(upper curve). Dashed line stand for a fit of a the experimental cell shape to Eq. 
(pS)) . which allow to obtain the parameters used in the numerical simulations. Dotted 
line is a translation of the dashed line to stress the self-consistency of the parameter 
extraction, b) Fit to Eq. of 7 valleys extracted from top curve in Fig. [T^. Red 
solid line stand for a dashed and dotted line in Fig. [T]i. c) Fit to Eq. ((29|) of 8 valleys 
extracted from bottom curve in Fig. [T]i. Blue solid line stand for a dashed and dotted 
line in Fig. [T^. 

task of determining the necessary and sufficient conditions for coarsening [16j, as well 
as the possible values (universality classes) for the coarsening exponent n, remains to 
be completed, lying outside the scope of this paper. Note that, once coarsening occurs, 
the fact that it interrupts [16] or not, and the value of exponent n, both depend on 
the type of leading and subleading nonlinearities that appear in the evolution equation. 
Hence, these provide additional criteria in order to propose a continuum description 
for a given physical system, although it can nevertheless be the case that two different 
subleading terms lead to the same coarsening exponent value, such as is the case for 
A2 and A210, both of them yielding n = 1/2 [331 EI]- Moreover, system dimensionality 
and the (possibly related) relevance of noise may influence the value of the coarsening 
exponent in the corresponding universality classes. 

In summary, we have provided a systematic method to obtain and classify the 
relevant nonlinear terms in interfacial systems which present morphological instabilities 
and the shift symmetry. Besides, we have applied our method to some reported systems 
in the literature as well as to novel experiments on ion plasma erosion performed for 
this purpose. 
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Figure 2. a) Ion plasma eroded Si(lOO) at <e = 3 h. Inset: Numerical simulation of 
Eq. ^ a.tts = 35 with v = I, K = 0.11, A210 = 0.052 and A = -0.05. The lateral 
size of the image is 50 /Ltm. b) Ion plasma eroded Si(lOO) at <e = 6 h. Inset: Simulation 
as in a) for — 70. The lateral size of the image is 50 /im. 
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